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Appendix A: 

This appendix contains the method developed in [5] for computing the 9-point DCT 
given by 

9-1 

A(n) = ^a(k) cos 

*=0 



2^9 



(2k + \)n 



for 0<n<8 



Method for the 9-point DCT: 

c x = -0.866025, c 2 =0.939693, c 3 =-0.173648, c 4 =-0.766044 

c s = 0.5, c 6 = -0.342020, c 7 = -0.984808, c 8 = -0.642788 

d x = a(3) + a(5), d 2 = a(3) - a(5), rf 3 = a(6) + a(2), tf 4 = a(6) - a(2) 

d s = a(l) + a(7), rf 6 = a(l) - a(7), rf 7 = a(8) + a(0), rf 8 = a(8) - a(0) 

</ 9 = a (4) + J 5 , d 10 = d x +d 3 ,d u =d xa +d 1 ,d x2 =d 3 -d 1 ,d u =d x -d 1 , d X4 = d x -d 3 

d x5 =d 2 -d 4 ,d x6 =d l5 +d s ,d xl =d 4 +d s ,d iS =d 2 -d !t ,d X9 =d 2 +d 4 

m x = c x d 6 ,m 2 = c 5 d 5 ,m 3 =c 5 d u ,m A = c 2 d x2 ,m 5 = c 3 d u 

m 6 =c 4 d X4 ,m 7 = c x d X6 ,m s = c 6 d xl ,m 9 = c 1 d n ,m w = c % d x9 

d 20 = a(4) - m 2 , d 2X = rf 20 + m 4 ,d 22 = <* 20 - w 4 , d 23 = d 20 + m 5 

d 24 =m x + m s ,d 25 =m x - m % ,d 2fl =m x +m g 

A 0 =d 9 + d n ,A x = m xo - d 26 ,A 2 =m 6 - d 2X ,A 3 = «„ 4« = d 21 - m 5 

A 5 =d 25 -m 9 ,A 6 =m 3 -d 9 ,A n = d 24 +m x0 ,A s =d 23 +m 6 

Thus, the 9-point DCT requires only 10 multiplications and 34 additions, respectively. 
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Appendix B: Source Code 

//*************************************** 
// File name: mdct.cpp 

#include <iostream,h> 
#include <ioraanip.h> 
#include <math.h> 
#include <time.h> 

#pragma hdrstop 
#include <condefs.h> 
#include "dct_lib.h M 

USEUNIT {"dct_lib, cpp") ; 

// 

int borg[100036] ; 
double 

bin [100036] , *ptbin f boutO [ 36] , boutl [36] ; 

// 

#pragma argsused 

int main (int argc, char* argv[]} 

{ 

int k,n; 
int i , cnt ; 
double diff; 
clock_t ckl,ck2; 

cout«" welcome to mdct test\n" ; 
build_costab ( ) ; 

// the random data 
randomize { ) ; 
for (k=0;k<100036;k++) 
borg[k] =random{255) ; 

for (k=0;k<100036;k++) bin [k] = (double) 

borg [k] ; 

// verify the MDCT and the fast methods 

cnt-100000; 

ptbin=bin; 

cout«"Verifymg MDCT3618 and MDCT3618F, 
cnt="«cnt<<endl<<endl ; 
f or (i=0; Kent; i + + ) 
{ 

MDCT3618 (ptbin, boutO) ; 
MDCT 3 61 8 F (ptbin, boutl ) ; 
ptbin++; 

for (n=0;n<18;n++) 
{ 

diff=( (boutO[n]- 
boutl[n] )<. 0000000001) 70:99999; 
if (diff>0) 

cout«diff<<setw(10) «bout0 [n] «setw (10) < 
<boutl [n] «"\n"; 
} 

} 

// verify the IMDCT and the fast methods 
ptbin=bin; 

cout«"Verifying IMDCT1836 and 



IMDCT1836F, cnt="«cnt«endl«endl; 
for (i=0; Kcnt;i++) 
{ 

IMDCT1836 (ptbin,bout0) ; 
IMDCTl836F{ptbin, boutl) ; 
ptbin++; 

for (n=0;n<36;n++) 
{ 

diff-( (bout0[n]- 
boutl [n] )<. 0000000001) ?0 : 99999; 
if (diff >0) 

cout«dif f «setw ( 10) «bout0 [n] «setw (10 ) «bo 
utl[n]«"\n"; 

} 

} 



cnt=100000; 

// computation time of MDCT 3 61 8 
cout« "Computation time of MDCT3618 
ptbin^bin; 
ckl=clock { ) ; 
for (i=0; Kent; 
{ 

MDCT3618 (ptbin, boutO) ; 
ptbin++; 

} 

ck2=clock() ; 

cout«setw ( 6) « ( ck2-ckl ) «endl ; 



// Computation time of MDCT3618F 
cout«" Computation time of MDCT3 618F 
ptbin=bm; 
ckl=clock ( ) ; 
for (i=0; Kent; i++) 
{ 

MDCT3618F (ptbin, boutl) ; 
ptbin++; 

} 

ck2=clock { ) ; 

cout«setw ( 6 ) « ( ck2-ckl ) «endl ; 



//Computation time of IMDCT1836 
cout«"Computation time of IMDCT1836 : "; 
ptbin=bin; 
ckl=clock() ; 
for (i=0; Kent; 
{ 

IMDCT1836 (ptbin, boutO) ; 
ptbin++; 

} 

ck2=clock () ; 

cout«setw(6}« (ck2-ckl) «endl; 



// Computation time of IMDCT1836F 
cout« n Computation time of IMDCT1836F : "; 
ptbin=bin; 
ckl=clock ( ) ; 
for (i=0; Kent; 
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{ 

IMDCT1836F (ptbin, boutl ) ; 
ptbin++; 

} 

ck2=clock ( ) ; 

cout«setw(6) « (ck2-ckl) «endl«endl; 

return 0; 
} 

//*************************************** 
// File name: dct_lib.cpp 



#include <iostream.h> 
#include <math.h> 
#pragma hdrstop 

#include "dct_lib.h" 

// 

#pragma package (smart_init) 

// some pre-computed consine tables 
double Cos99[9] [9] ; 

double Cos99W[9]; // used for DCT99W 

double Cosl818[18] [18]; 

double Cosl818F[18] ;// DCT1818F 

double Cos3618 [18] [36] ; // shared for 

mdct3618, imdctl836 

double Cos3618F[18] ; 

double Cosl836F[18] ; 



//===ra======m =====™===============:- 

// the COS tables for various size 

void build_costab (void) 
{ 

int n,k; 

// DCT99 

for (n=0;n<9;n++) 

for (k=0;k<9;k++) 

Cos99 [n] [k]=cos (M_PI/18* (2*k+l) *n) ; 

// for DCT99W 
Cos99W[l]=Cos99[l] [7] ; 
Cos99W[2]=Cos99[2] [0] ; 
Cos99W[3]-Cos99[2] [2] ; 
Cos99W[4]=Cos99[2] [3] ; 
Cos99W[5]=Cos99[2] [1]; 
Cos99W[6]=Cos99[l] [5] ; 
Cos99W[7]=Cos99[l] [8] ; 
Cos99W[8]-Cos99[l] [6] ; 



// DCT1818 

for (n=0;n<18;n++) 

for (k=0; k<18; k++) 

Cosl818 [n] [k]=cos {M_PI/3 6* (2*k+l) *n) 
// DCT1818F 



for (k=0;k<18;k++) 

Cosl818F[k]=2*Cosl818[l] [k] ; 

// DCT3618 

for (n=0;n<18;n++) 

for (k=0;k<36;k++) 

Cos3618[n] [k] =cos (M_PI /72* ( 2*k+19 ) * (2*n 

+D); 

// DCT3618F 

for (k=0;k<18;k++) 

Cos3618F[k]=2*cos (M_PI/72* (2*k+l) ) ; 

// IMDCT1836F 
for (k=0;k<18;k++) 

Cosl836F[k] =cos (M_PI/72*19* <2*k+l) ) ; 

} 



// 9-9 DCT / standard 

void DCT99 (double *a, double *A) 
{ 

int n,k; 
double s; 

for (n=0;n<9;n++) 
{ 

s=0; 

for (k=0;k<9;k++) 
{ 

s+=a[k] *Cos99 [n] [k] ; 

} 

A[n]=s; 
} 

} 



// 9-9 DCT, Winogard 

= = 

void DCT99W (double *a, double *A) 

{ 

double d[64] ,m[64] ; 

d[l]=a[3)+a[5] ; 
d[2]=a[3]-a[5]; 
d[3]=a[6]+a[2] ; 
d[4]=a[6]-a[2] ; 
d[5]«a[l]+a[7] ; 
d[6]=a[l]-a[7] ; 
d[7]=a[8]+a[0] ; 
d[8]=a[8] -a[0] ; 

d[ 9]=a[ 4]+d[5J; 
d[10]=d[ l]+d[3]; 
d[ll]=d[10]+d[7]; 
d[12]=d[ 3]-d[7]; 
d[13]=d[ l]-d[7]; 
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d[14]=d[ l]-d[3]; 
d[15]=d[ 2]-d[43; 
d[16]=d[15]+d[8] ; 
d[17]=d[ 4]+d[8]; 
d[18]-d[ 2]-d[8]; 
d[19]=d[ 2]+d[4]; 

m[ l]=Cos99W[l]*d[ 6] 
m[ 2)=Cos99W[5] *d[ 5] 
m[ 3]=Cos99W[5]*d[ll] 
m[ 4]=Cos99W[2j*d[12] 
m[ 5]=Cos99W[3] *d[13] 
m[ 6]=Cos99W[4]*d[14] 
m[ 7]=Cos99W[l] *d[16] 
m[ 8]=Cos99W[6]*d[17] 
m[ 9]=Cos99W[7]*d[18] 
m[10]=Cos99W[8]*d[19] 

d[20]=a[ 4]-m[2]; 
d[21]=d[20]+ra[4] ; 
d[22]=d[20]-m[4] ; 
d[23]=d[20]+m[5] ; 
d[24]=m[ l]+m[8]; 
d[25]=m[ l]-m[8]; 
d[26]=m[ l]+m[9]; 

A[0]=d[ 9]+d[ll]; 
A[lJ«m[10]-d[26] ; 
A[2]=m[ 6]-d[21]; 
A[3]=m[ 7]; 
A[4]=d[22]-m[ 5] ; 
A[5]=d[25]-m[ 9] ; 
A[6]=m[ 3]-d[ 9] ; 
A[7]=d[24]+m[10] ; 
A[8]=d[23]+m[ 6] ; 
} 



// 18-18 DCT, standard 

/ / ^ == =™==============™==-™== 

void DCT1818 (double *a, double *A) 
{ 

int n,k; 
double s; 

for (n=0;n<18;n++) 
{ 

s=0; 

for (k=0;k<18;k++) 
{ 

s+=a[k]*Cosl818[n] [k] ; 

} 

A[n] =s; 
} 

} 



//================= 

// 18-18 DCT, fast 



void DCT1818F (double *x, double *X) 
{ 

int n,k; 

double a [9] ,b[9] ,A[9] ,B[9] ,Bp[9] ; 

for(k=0;k<9;k++) 
{ 

a[k]=x[k]+x[18-l-k] ; 
b[k]=(x[k]-x[18-l-k] ) *Cosl818F[k] ; 
} 

DCT99W (a, A) ; 
DCT99W(b / Bp) ; 

B[0}=Bp[0]/2; 

for(n=l;n<9;n++) B[n]=Bp[n] -B[n-1] ; 

for (n=0;n<9;n++) 
{ 

X[2*n]=A[n] ; 
X[2*n+l]=B[n] ; 
} 

} 



// 36-18 DCT, standard 
//===================™===—======= 

void MDCT3618 (double *a, double *A) 
{ 

int n,k; 
double s; 

for (n=0;n<18;n++) 
{ 

s=0; 

for (k=0;k<36;k++) 
{ 

s+=a[k]*Cos3618[n] [k] ; 

} 

A[n]=s; 
} 

} 



// 36-18 DCT, fast 

//= _™ ======:= _ == ^== = == = =™============== 

void MDCT3 618F (double *y, double *Y) 
{ 

int n,k; 

double x[18] , Yp[18] ; 

for (k=0; k<9;k++) 

x[k]=(-y[26-k] -y[27+k] ) *Cos3618F[k] ; 
for (k-9;k<18;k++) 

x[k]=( y[k-9 ]-y[26-k])*Cos3618F[k]; 

DCTl818F(x, Yp) ; 
Y[0]=Yp[0] /2; 

for (n=l;n<18;n++) Y [n] =Yp [n] -Y [n-1 ] ; 
} 
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File name : dct_l ib . h 

//*************************************** 



// 18-36 IMDCT, standard 

//========== ============================ 

void IMDCT1836 (double *Y, double *y) 
{ 

int n, k; 
double s; 

for (n=0;n<36;n++) 
{ 

s=0; 

for <k=0;k<18;k++) 
{ 

s+=Y[k] *Cos3618 [k] [n] ;//share 
with the same table as 3618 
} 

y [n]=s; 
} 

} 

//=== _ =TO= ^ ======= „==================== 



#ifndef lib_dctH 
#define lib_dctH 

//==================================== 

void build_costab (void) ; 
void DCT99 (double *a, double *A) ; 
void DCT99W {double *a, double *A) ; 
void DCT1818 (double *a, double *A) ; 
void DCT1818F (double *a, double *A) ; 
void MDCT3618 (double *a, double *A) ; 
void MDCT3618F (double *a, double *A) ; 
void IMDCT1836 (double *Y, double *y) ; 
void IMDCT1836F (double *Y, double *y) ; 
#endif 



// ===== ==^=^==================-====== == 

// 36-18 DCT, fast 

// == =======™=============-============— 

void IMDCT1836F (double *Y, double *y) 
{ 

int n, k; 

double Yp[18] , yppp [18] , yp[36] , s; 

for (k=0; k<18;k++) 

Yp[k]=Y[k]*Cos3618F[k} ; //the same 

table 



DCTl818F(Yp, yppp) ; 

for(n= 0;n< 9;n++) 
for(n= 9;n<10;n++) 
for(n=10;n<27;n++) 
for(n=27;n<36;n++) 



yp [n] = yppp [n+9] ; 
yp[n]= 0; 

yp[n]=-yppp[27-n] ; 
yp[n] =-yppp[n-27] ; 



s=0; 

//for(k=0;k<18;k++) s+=Yp[k]; 
for(k=0;k<18;k++) s+=Y[k] *Cosl836F[k] ; 
y[0]=s; 

for (n=l;n<36;n++) y [n]=yp [n] -y [n-l] ; 
} 
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